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Abstract 

We investigate the possibility of generalizing differential renormalization of 
D.Z.Freedman, K.Johnson and J.I.Latorre in an invariant fashion to theories with 
infrared divergencies via an infrared R operation. Two-dimensional a models and 
the four-dimensional </> 4 theory diagrams with exceptional momenta are used as 
examples, while dimensional renormalization serves as a test scheme for com- 
parison. We write the basic differential identities of the method simultaneously 
in co-ordinate and momentum space, introducing two scales which remove ul- 
traviolet and infrared singularities. The consistent set of Fourier-transformation 
formulae is derived. However, the values for tadpole-type Feynman integrals in 
higher orders of perturbation theory prove to be ambiguous, depending on the 
order of evaluation of the subgraphs. In two dimensions, even earlier than this 
ambiguity manifests itself, renormalization-group calculations based on infrared 
extension of differential renormalization lead to incorrect results. We conclude 
that the extended differential renormalization procedure does not perform the 
infrared R operation in a self-consistent way, as the original recipe does the 
ultraviolet R operation. 
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1 Introduction 



Differential renormalization |IJ can successfully perform an ultraviolet R operation 
by means of replacing essentially singular quantities in co-ordinate space with deriva- 
tives of less singular expressions, and after that, by integrating the derivatives by parts 
in order to turn them into powers of external momenta. Then, at first sight, one sees 
no serious obstacles — it would just seem very natural indeed — that the same trick 
with singular expressions in momentum space could likewise automatically perform an 
infrared R operation 0. 

A need of the latter in renormalization-group calculations (leaving alone operator- 
product and other asymptotic expansions) may arise of two causes. In the first place, 
a quantum field theory may be intrinsically infrared- dangerous owing to the fact that 
the theory with zero masses may not strictly exist at all as, for example, a massless 
scalar field in two dimensions. However, real computations with nonzero masses are 
always exceedingly cumbersome technically while, on the other hand, we know that 
in the so-called "massless" renormalization schemes || (like the minimal subtraction 
scheme |4], || of dimensional or analytic regularization) there is no essential dependence 
of renormalization-group functions on dimensional parameters (except for a possible 
mixing in renormalizations of the masses themselves). Thus, we would like to set them 
all to zero from the very beginning but use an infrared R operation to avoid the false 
infrared singularities. 

In the second case, artificial infrared divergencies may be caused by a singular choice 
of external momenta of a Green's function, quite well defined in a general situation. For 
technical reasons, such a choice may seem more convenient, or it may even be the only 
available choice that would allow us to compute a complicated multiloop vertex-type 
graph analytically by means of an infrared rearrangement ||. 

Although there are no theoretical grounds, like the causality principle of the local 
quantum field theory, for attributing a physical meaning to the infrared R operation 
in itself — the reasons look quite technical both with intrinsic and artificial infrared 
divergencies — however, in many practically important cases its use may be crucial for 
achieving computational results [7L H. 

Introducing the infrared R operation into the formalism of differential renormaliza- 
tion, we should keep track of three important points. First of all, the scheme should 
retain its self-consistency: the results should never depend on the order of manip- 
ulations we perform. The second requirement is that we should distinguish two di- 
mensional scales: a new infrared renormalization scale is, and an ultraviolet scale \i 
already present in the formalism. Otherwise, infrared logarithms would be the same as 
ultraviolet logarithms, and we would never be able to extract from the finite renormal- 
ized Feynman diagrams the correct information about renormalization-group functions, 
contained in their /i dependence. And last but by no means least point is the invariance 
of the regularization ||, which ensures the fact that formal symmetries of a quantum 
field theory are retained after renormalization || [Uj] . 

We formulate a natural generalization of differential renormalization to theories 
where, besides usual ultraviolet singularities, intrinsic or artificial infrared divergencies 
do appear. Our particular attention is paid to two-dimensional a models and to the 



four-dimensional 4 theory. Throughout our renormalization-group calculations, we 
compare the results of the extended differential renormalization with the minimal sub- 
traction scheme of dimensional regularization as an example of an undoubtedly valid, 
though different, renormalization scheme in the range of scheme arbitrariness. 



2 D=2: Intrinsic infrared divergencies 



As the first example, we consider two-dimensional a models [IT], 0] where we meet 
infrared difficulties of both intrinsic and artificial origin. 

The starting point of differential renormalization in a two-dimensional theory is the 
following replacement which effectively subtracts an infinite local ultraviolet counter- 
term proportional to the 5 function in co-ordinate space: 



1/x 2 
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where D x denotes (d/dx^) 2 (we work in Euclidean space). The other way of looking at 
the matter is to consider eq. (|TJ) as a real equality under a regularization, bearing in mind 
that the regularized expressions are renormalized afterwards by minimal subtractions 
which simply turn the regulator M— >oo into a finite renormalization scale. 

With the aid of eq.(|IJ), integrating the derivatives by parts and abandoning surface 
terms (which is one more prescription of the regularization), we can derive a basic set of 
formulae for Fourier transformations. We start with the ultraviolet-divergent integrals 
of the form 

ln n (a; 2 M 2 ) 

-■ (2) 
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Replacing 1/a; 2 through eq. ([[]), we integrate one d/dx^ by parts. Of course, we suppose 
the Leibniz differentiation rule to remain valid under the regularization. The first 
derivative of the logarithm is a nonsingular function in two dimensions, and hence, it 
can be substituted explicitly as 2x^/x 2 without any possible 5-function-like terms. The 
resulting equation can be solved with respect to I n : 
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The integral on the right-hand side has no dangerous singularities, and we can unam- 
biguously evaluate it, using an intermediate analytic regularization which is taken off 
in the end, 
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F m are mth derivatives of that function at zero, 



F = l, F 1 = F 2 = 0, F 3 = 4C(3), F 4 = 0, F 5 = 48 C(5), 
F 6 = 160C 2 (3), F 7 = 1440 C(7), 

C(m) is the Riemann zeta function, and C is the Euler constant. 

Since the initial integrals I n , eq.(Q), were logarithmically divergent, in the most 
general situation we could add to the renormalized expressions (f|) some constants 
C n , reflecting a local arbitrariness. In the context of differential renormalization, we 
could leave an explicit 5-function term in eq. ([[]). That would simply be equivalent 
to redefining the renormalization scale M. Therefore, we have chosen to absorb the 
arbitrariness into the definition of M, and no more C n could ever appear with the 
minimal subtractions on eq.(§). 

Now, we can try to execute precisely the same procedure for infrared rather than ul- 
traviolet divergencies, starting from momentum space with an infrared regulator mass 
v— >0, which becomes another finite renormalization scale after the minimal subtrac- 



tions: 



l/p 2 = §□„ InV/z, 2 ) , (6) 
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d 2 p e-*» % In V/, 2 ) = -2vri % (-)« ± ( " ) F n . m ln m (x 2 iV 2 ) , (7) 

P n + 1 m=Q \ m J 

Again, we have denoted N 2 = ~ v 2 e 2C . (Everywhere, M 2 and N 2 accompany x 2 , while 
p 2 is divided by jj? or z/ 2 ). 



2.1 The simplest example 

The fact that there really is an infrared difficulty for differential renormalization in 
two dimensions is quite obvious from considering the simplest one-loop tadpole diagram 
which appears in a models (analyzed below in subsects.|2T4] and |2~5|) : 



(9) 



The original method of ref.ffl would simply set this diagram to zero, leaving no trace 
of any ultraviolet logarithm. Therefore, in order to find a correct contribution to the 
(3 function, we need a special procedure for defining such integrals which involve no 
external momenta or co-ordinates and are both ultraviolet- and infrared-divergent. 

Let us start with the outline of our ideas in application to this integral as a simple 
example. We always assume that the whole scheme can be treated as an invariant 
regularization procedure which allows formal manipulations with divergent integrals. 
Then, we can identically transform eq.([|) to introduce an auxiliary external momentum 
argument k 
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The first integral on the right-hand side of eq.flTOD is infrared- divergent but ultraviolet- 
regular, while the second integral has only an ultraviolet singularity. Now we use the 
Fourier transformation to rewrite eq.flTUP as 



d2x ^e ikx (k Jd 2 p e -^- - U 2 pe-^^) fd 2 qe-^%. (11) 
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One Fourier integral is regular [it is the simplest particular case of eq.([7|)] 



/d 2 pe-^^ =-2vri^. (12) 
J v x 
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The other Fourier integral in eq.flTTD requires an infrared renormalization. It is io, 
eq-(i, 

e -ipx 



d 2 p^r- = -7T ln(x 2 iV 2 ) . (13) 
P 

In the extended differential renormalization, the latter equation is derived from the 
momentum-space differential identity (§) via integration by parts and abandoning sur- 
face terms which are absorbed into local p-space counterterms in the formalism of the 
infrared R operation . The result involves a new infrared renormalization scale. For 



the case of the Fourier transform of a free scalar-field propagator (13) , the physical 
meaning of the R operation can easily be understood: it recovers the most essential 
(logarithmic) part of the massive propagator that remains as the mass goes to zero 

e~ ipx 

d 2 p -s n = 27T Kq {m\x\) = — 7r In ( \ m 2 x 2 ) + o(l) , m — > . 

p + m v * / 

In this interpretation, the infrared differential regularization parameter in eq.(^) is pro- 
portional to the mass v 2 =m 2 e~ 2C . However, the general R formalism is not reduced to 
introducing a mass. It renormalizes any infrared- divergent diagram in a dual way to the 
ordinary R operation with local x-space counterterms, the structure being completely 
analogous, but with the x^*p interchange. 



To finish the computation of the tadpole diagram (|9|), we substitute eqs.(|T2|) and 
(|13|) into eq . ([TT]) . The integrals over x are taken by the original rules of differential 
renormalization, which lead to eqs.(|3]) and (f|), in particular, 

r -r h h 2 r p ikx h 2 

/ d 2 x e ife - % ln(a; 2 iV 2 ) = -2th % In ^ , / d 2 x ^- = -tt In ^ , 

■/ x k v J x /i 

where the ultraviolet scale [i has appeared. In the resulting expression for the diagram 
(H) the dependence on the auxiliary momentum cancels, and we obtain n ln(/i 2 / u 2 ) . 
Differentiating this with respect to ln/z 2 , we reproduce the correct contribution to the 
(3 function. 

A systematic study of the implications of the infrared differential renormalization 
is presented below. First, in subsect. |2~2| , we try to derive a maximum set of relations 
(including singular functions) from the requirements of invariance of regularization. 
Since the maximum set entails a contradiction (subsect{2~3|), we restrict ourselves to 
avoid any explicit singular differetiation. However, even without that, we recover all 
the essential consequences, so that the more cautious approach turns to be equivalent 
to the unrestricted one. 

2.2 The consistent set of Fourier integrals 

Having added the momentum-space differential identity (H), we ought to check the 
consistency of the scheme as a whole. We expect of a good regularization that it is 
invariant |J, that is, permits us to do with our integrals some formal manipulations, 
which is important for retaining symmetries of quantum field models. One of those 
features is the possibility of differentiating the integrals with respect to external pa- 
rameters, for example, the Fourier-transformed argument, to compute other integrals, 
with external Lorentz indices, as we have done in deriving eq.(|3|). Other natural re- 
quirements are the possibilities of changing the order of integrations, of shifting integra- 
tion variables, and of identically transforming subintegral expressions, like expanding 
(p—q) 2 or canceling the numerator with the denominator. Such transformations are 
often performed, in particular, when one does superfield algebra in supersymmetric 
models to reduce supergraphs to usual Feynman integrals. 

Believing that the extended differential regularization can be treated as an invariant 
overall regularization, let us apply —O x to i n to obtain, after canceling the denominator 
p 2 in eq.(H), new integrals 

j n = -D x i n = j d 2 p e-** \n n (p 2 /u 2 ) , (14) 

which contain no infrared divergencies, as differs from eq.(^). Note also that any 
arbitrary renormalization-specific constants c n added to i n are annihilated by differen- 
tiation, so that j n will in any case be determined uniquely. 

At ?2=0, we must evidently get jo=47r 2 82(x). On the other hand, differentiating the 
right-hand side of eq.(|8]) for io, produces the laplacian of the logarithm. Hence, 



O x ln(x 2 N 2 ) = 4vr 8 2 {x) 



(15) 



Of course, in fact, eq.(|I~5D is true without any reference to differential regularization (N 
drops out, too). The identity can as well be derived by some other means, for example, 
by explicitly introducing the causal rule of passing over x 2 =0, expanding the left-hand 
side of eq. (|15D , and directly evaluating the integral of it. The extended differential 
regularization reproduces the result correctly. 

Now we proceed to n—1. The laplacian applied to the logarithm squared in ii, 
eq.(||), gives us two terms: when both derivatives act on one logarithm, we use eq. ([l5|) , 
while the nonsingular first derivatives of the logarithm can be substituted directly. 
Thus, we get 

ji = -D x i x = -An 2 6 2 {x) ln(x 2 N 2 ) -An/x 2 . (16) 



To attach a meaning to eq . (jig ) , we need to extend the definition of the logarithm to 
zero point in co-ordinate space. However, our set of formulae already includes such a 
definition: if we expand, on the right-hand side of eq.(|l|), the action of the laplacian, we 
get besides 1/x 2 an additional term which equals nS 2 (x) In (x 2 M 2 ) in view of eq. (TT5|) . 
We have decided that by definition of M this local term is absent in eq. ([!]); therefore, 



5 2 {x) \nx 2 = 6 2 {x) ln(l/M 2 
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It is this formula that shows us explicitly what the regularization really changes. For 
the rest, everything looks as if nothing has been regularized at all, provided we manage 
to take proper care of singular functions. 

As soon as we have established eqs. ([T7|) and (|T5|) , we can uniformly obtain all the 
integrals j n , eq.(IIl), from i n , eq.(W), arriving at the general formula 
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Quite analogously, acting with — O p on eq.(Q) and remembering eq.(|^), we deduce 
the x<-^p complementary formulae 
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Evidently, the ratios N 2 /M 2 in eq.([l8D and v 2 / 'fi 2 in eq.([20|) are the same. 

Another observation is that, as J n is ultraviolet regular, it does not generate an 
ultraviolet scale parameter by itself. Therefore, M in eq.(pO[) is just a dummy argument 



(strictly related to /i on the right-hand side) while the real stamp of the regularization 
on J n is the infrared scale v from eq . (|19D . The self-consistency of replacing M in 
eq.tJZUp can be checked directly by re-expanding the logarithm with an arbitrary mass 
argument A by the Newton's binomial, and resumming the right-hand side. The result 
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actually looks the same as eq. fl20D . Likewise, we can freely change v in eq.(|l8|) to a 
dummy variable A which can even be set equal to /x, 



J d 2 p e- [px ln n (p 2 /A 2 ) = 4tt 
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The intermediate integrals with Lorentz indices, eqs.([3D and (0), are both ultraviolet- 
and infrared-regular. Thus, they introduce no scales at all, \i and v being their dummy 
arguments. 

The next step in checking the consistency is to take the Fourier transforms of 
eqs. (|2l"l) and (|22|) . One easily sees that the resulting formulae allow us to recursively 
recover the very first integrals: /„, eq.(0), and i n , eq.(^). The double circle closes. 
By induction, using binomial re-expansions, we can prove that formulae (Q) and (§) 
are exactly reproduced without any C n and c n , which reflects the fact that we made 
eqs.(Ti7|) and ([19]) agree with the differential relations (|1|) and (|). For other masses 
under the logarithms, we get quite definite additive terms, 
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with A 2 = 4 A 2 e" 2C in eqs.(P)-P 



2.3 Ambiguities in repeated integrals 



Everything looks fine at the moment, and we could start to compute diagrams. 
However, the point is that we actually need a little more than just switching to the 
Fourier transforms and back. Namely, we expect that a property similar to associativity 
of convolution will hold. Suppose, there are three quantities 



a n (p) = J d 2 x e ipx A n (x) , A n {x) 
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(ra=l,2,3), and we are going to evaluate the integral of their product in momentum 
space. We can represent it in three different ways, substituting eq. (p5|) either for a 1; 
a 2 , or a 3 : 



X = I d 2 p ai(p) a 2 (p) a 3 (p) 

d 2 x Ai(x) J d 2 p e~ ipx a 2 (p) a 3 (p) 
d 2 x A 2 (x) J d 2 p e~ ipx ai(p) a 3 (p) 
d 2 x A 3 (x) J d 2 p e- [px Gl (p) a 2 (p) , 
where the integral over x can be evaluated by a trick of inserting a unit 
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and subsequently applying eq.([L9]). Either way (|26|)-(p8|) should lead to the same result 
in the framework of a consistent regularization. Let us, however, substitute particular 
expressions 
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All the integrals that arise can be evaluated with the aid of eqs.(^l])-(|24]), after products 
of different logarithms are re-expanded in powers of one of them. The results are 
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We see that there are two ways of computing the integral, which give us two different 
answers. The resulting expressions (|3~ll) and (|32]) coincide only if we set z/ 2 =/x 2 , that is, 
use the same scale for renormalizing both ultraviolet and infrared divergencies. As has 
been said at the beginning, this would not be satisfactory for renormalization-group 
applications. 

One may suspect that the inconsistency has resulted from undue liberty in our 
dealing with singular functions. Indeed, in establishing our set of formulae, we exceeded 



the necessary minimum of operations, essential for evaluating Feynman integrals. So, 
in principle, we can do without explicit differentiation the result of which might involve 
a singular function. 

Then, we ought to exclude all the singular relations, fljjP, (|T6|), (|17D, and (|19l), thus, 
being unable to get j n , eq. flllD , from i n . Nevertheless, the basic divergent integrals 
themselves I n and i n , eqs.(£|) and (|[), can be derived from the differential identities 
eqs.(|l|) and (||) through intermediate finite integrals, eqs.(0) and as before. Re- 
expanding the logarithms, we reproduce the integrals with arbitrary mass arguments, 
); and finally, eqs. (p2|) and (|2TD are obtained as their resummed Fourier 
transforms with the particular cases of eqs.(|18|) and (|20|). 

Although the trick of eq. (|29|) becomes of no use, we can, all the same, define the 
general external-parameter-independent integrals in eqs.(PB|)-(P^) by utilizing the as- 
sumed invariance of the regularization: 
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where we have made use of eqs.(||), (^) and (0). Note that the term with m=0 is 
always present, even at A 2 =iV 2 . In the same way, we can evaluate the corresponding 
momentum integrals 
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Now, for the chosen ai, a 2 , and 03, as in eq.(|30|), the integrals on the right-hand 
sides of eqs.(f26|)-(|28D are directly calculable, and we reproduce the two different results 
(HD and © again. 



2.4 N=2 supersymmetric sigma model 

Yet, let us see how grave this contradiction is from the practical viewpoint of 
computing Feynman diagrams. We consider first the four-loop approximation in the 
N=2 supersymmetric two-dimensional a model. After the supergraph algebra in the 
background-field formalism is done [|EJ], there appears just one nontrivial divergent 
four-loop momentum integral which is presented in fig.[l]a. The rest of the divergent 
integrals are of the primitive tadpole type (fig.|l|b); we leave them alone for the moment. 




(a) (b) 

Figure 1: The divergent four-loop momentum integrals in the N=2 supersymmetric 
two-dimensional a model. Lines represent 1/p 2 propagators, arrows stand for additional 
Pn in the numerator, and small dots denote contractions of Lorentz indices. 



The model involves intrinsic infrared divergencies due to zero masses on the lines 
without components of momenta in numerators. Besides, there may appear artificial 
infrared divergencies if some of the external momenta in fig.|l|a are set to zero. In 
fact, for technical reasons, we simply cannot do without this, to reduce the problem 
to successively evaluating one-loop propagator-type graphs. The three external lines 
of the diagram correspond to identical operator structures of the background field. 
Therefore, nullifying one of the three momenta, we obtain two different graphs shown 
in fig.0, while setting all of them to zero, we get a pure tadpole. The diagram of fig.^b 
and the complicated tadpole graph contain artificial infrared divergencies. 




w 00 

Figure 2: Two possibilities that arise when one of the external momenta is set to zero. 



We can easily compute the two one-loop subgraphs both in co-ordinate and in 
momentum space fl25|) by means of eqs.(0), (|^), ([ 



and (HI): 
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Now, the diagram of fig.^a can be calculated as follows. We first compute its two-loop 
propagator-type subgraph in momentum space as the square of the one-loop expression 
Then, we transform the subgraph to co-ordinate space 
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The whole diagram is just a product of eq.(]3"7|) and the infrared one-loop subgraph 
(|36|) in x space. Re-expanding the infrared logarithm squared in terms of ln(x 2 M 2 ), 
and switching back to momentum space by means of eq.(HI), we arrive at the result 
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The same result is obtained if we leave the one-loop subgraph as it is, but reduce 
eq.(|3"TD to ln(x 2 iV 2 ), take the Fourier transform of the product via eq. fl2"B|) , and finally, 
re-expand it in the ultraviolet logarithms. 

In this way, using only the Fourier transformation formulae, we can compute any 
diagram of recursively one-loop propagator structure. In particular, we consecutively 
evaluate 
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The answer for the tadpole diagram at zero momenta in figfja proves to be ambigu- 
ous. If we first compute its infrared-regular two- loop subgraph (|3"7|), and then evaluate 
the tadpole by the ^-function trick ( p9|) or by eq. (|33|) , we get 




(4tt)< 



±lnV//^)+fC(3) Hu 2 /^ 2 ) 



(41) 



like in eq.(|3T]). The use of the direct momentum-space formula (g4|) leads to the same 
answer. On the other hand, starting from the infrared-divergent two-loop subgraph 



|H), we repeat eq.fl32j): 




(AnY 



ilnV/V)- |C(3) ln(^ 2 /^ 2 ) 



(42) 



Both answers flilD and ( f42"l) can be obtained from eqs.(^) and (0), respectively, by 
replacing p 2 with v 2 . 

To estimate the meaning of the results, it is worth comparing them with the un- 
conditionally consistent scheme of dimensional renormalization. The infrared R and 
ultraviolet R operation ought to be done explicitly then. Without entering into details, 
let us present as an example (fig.0) the structure of renormalizations for the diagram 
of fig.|b. 
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+ 2 
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Figure 3: The structure of the incomplete ultraviolet R' operation and the full infrared 
R operation. A subtracted ultraviolet subgraph is shrunk to a point, while an infrared 
subgraph is just deleted. After the R, we have omitted zero tadpole graphs. 



The subtraction operators K and K pick out poles in e=(2—D)/2 in momentum 
and co-ordinate space, respectively. We keep different ultraviolet and infrared renor- 
malization scales. Thus, every momentum loop integration is, as usual, accompanied 
by a (fi 2 ) £ factor, while the proper dimension of infrared counterterms proportional 
to 62-2e(p) is restored by additional denominators of (^ 2 ) e . For example, in fig.|3] we 
should substitute 

one (/x 2 ) £ coming from the integration annihilated by the S function. 

Instead of doing the infrared R operation in the diagram language, as in fig.0, it 
may technically be easier to subtract singular momentum-space expressions directly by 
means of a simple formula || 

-(n+l)e 



1 



R 



1+ne 



(?) 



1+ne 



+ 7T 



n 



leaving only pole terms on subtracting complicated subgraphs with a prefactor. 
Here are the results of dimensional renormalization in the so-called G-scheme 
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InV/V)- flnV/V) ln(^ 2 )+lnV/r) biV/V 
+ 8C(3) ln(p 2 //i 2 ) - f C(3) ln(z/ 2 //i 2 ) + |C(4) , 



(43) 



(47T) 
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+ | InV/M 2 ) - | ln 3 (p 2 /^) In^/V) + | ln(p 2 //i 2 ) InV/V) 



| InV/ 

-|C(3)[ln(//^)+ln(r/7^) + § C(4) , 



(44) 



(4vr) 4 M 
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+ 



C(3) 



+ iln 4 (^ 2 )-4C(3) h^V/i 2 ). (45) 



Comparing eqs. (|38|)-(f42|) with eqs.(43)-(45), we come to the following conclusions. 
The results of the extended differential renormalization disagree with the dimensionally 
renormalized diagrams obtained by subtracting pole contributions, RR—(1—K)R'R. 
Moreover, this disagreement does not fit into the framework of the scheme arbitrariness 
by just redefining /i and v: eq.([P]) involves a ln(j? 2 /// 2 ) term absent in eq. (fTO|) . Forming 
the weighed sum, according to fig.0, does not help in improving the agreement. The 
tadpole graphs are not well-defined at all, none of the results coinciding with eq. (H5|) , 
neither eq.(|l|) nor (|42|). 

One cannot completely avoid tadpole-graph ambiguities even at [i=v (this would 
be sufficient for asymptotic expansions). Already for the five-loop tadpole diagram 
shown in fig.|], we can get three different answers at [i=v^ starting from its one-, two-, 
or three-loop infrared- divergent subgraphs. This is not just an artificial example: the 
graph appears in the five-loop calculations of ref.||. 



(4tt) 
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C(5), t-C(5). 



Figure 4: The five-loop tadpole graph, the value of which in the extended differential 
renormalization is ambiguous even at \i=v. 



Trying to extract from eqs. (j38|)-([42|) the contributions to the renormalization-group 
(3 function by differentiating the expressions with respect to In /i 2 , we find that the result 
depends on ln(i/ 2 //z 2 ) and ln(p 2 //i 2 ). In the present case, the only nontrivial diagram 
of fig.|l|a is proportional to a new invariant structure in the Riemann tensors and the 
background field, which was not present in lower orders of perturbation theory. Thus, 
no compensation of higher-order logarithms from lower-order graphs can occur, as it 
happens in ordinary renormalizable field theories. The ignored contributions from the 
primitive tadpole-reducible graphs of fig.|l|b cannot save the situation with higher-order 
logarithms of the momentum since they contain nothing but ln 4 (i^ 2 / /i 2 ). This problem 
is present in dimensional renormalization as well, if we attempt to find the (3 function 
by using the finite logarithms instead of applying the traditional method M based 



on the first-order pole in e. Even at p 



-v 



-fj, , fig.Q with the finite parts of eqs.( 



and (53) does not reproduce the correct (3 function |T2|] extracted from the pole term, 
although the dimensionally renormalized tadpole graph gives the consistent answer in 
both ways. The cause of these difficulties is non-renormalizability of the model, as a 



consequence of which the renormalization-group-like equations can only be written in 
terms of the pole contributions rather than directly in finite logarithms [14, 11 1 . 



2.5 A renormalizable sigma model 



It would be natural then to consider a renormalizable model, where the renormal- 
ization-group functions can as well be computed from finite logarithms. As such a 
model, we can choose the n-component n field in two dimensions, 



C = {2h)- 1 (^n 



11" 



which can be rewritten as a special case of the simplest bosonic a model |TT]] for n—1 
independent components 



h 



9jk - 5 jk + 1 _ j, 2 (j>j4>k 



For this particular form of the metric with the co-ordinates <pj, we have 

g jk = 5 Jh - h <\P<\>\ r jk l = h 9jk <f> 1 , 

Rjklm = h (gjl gkra ~ gjm gkl) , Rjk — [Tl — 2) h gjk ■ 

The invariant charge h=Z~ l h is determined by the charge renormalization constant Z. 



The contributing Feynman diagrams in the background- field formalism [11 up to two 
loops are presented in fig.||. 



-R jk (d^) d^ k 
+ i (2 R a 3 R ak + 3 R abc j R abck ) (d^) d^ k 



(b) 



\Rj a \ R ab (0^)8, 



(d) 



Figure 5: The one- and two-loop corrections to the effective action of the two- 
dimensional bosonic a model without torsion. The contributions to Z are obtained 
by normalizing to the tree term | (c^0 J ) gjk d^<j) k . 

The results of dimensional renormalization for the n field are 



R'R (fig.5a) = -(n - 2) g jh (d^)(d^ k ) 
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h m — 

e v 



(46) 



R'R (fig.|b) = i ( n + l)(n - 2) (d^)(d^ k 
R'R (fig.|c) = |(n-2)V(^')(^0 fc ) 



(47T) 
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f2'E (fig.5d) 



P 



{n-2)g jk (d^ j )(d v cf> k ) 



h 2 



(4tt) 



2 



+ ln z ^ - 2 In - ln z ^ - 2 In ^ + 1 



7/ 



P 2 



1 1 

? " e 

"2 



In 2 ^ - 2 In ^ + 2 



1/ 



(49) 
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These formulae lead to the renormalized 
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Independence of the invariant charge on fi 2 leads then to the following expression for 
the (3 function 



(3&\m{h) = h 



dZ 
d In /i 2 



1 -h—\ Z 

dh 



-1 



, h2 ( h \ 
- 2 < B - 2 >^ " (51) 



The result is the same as obtained from the 1/e-pole terms in eqs.(HH)-(^0). It also 
agrees with the known (3 function for the ri-field 0(ri) sigma model |1| where no 
infrared R operation was ever used. 

To be capable of evaluating the Feynman integrals of fig.[|(d,e) in the extended 
differential renormalization, we need to extend the set of formulae to a more general 
kind of expressions. The new differential relation is 



ln n (a; 2 M^ 



-in! — 

2 dx. 



1 » ln m (x 2 M 2 ) ^ , 

— 2^ ; +7r5 2 (x) 



X 



m=0 



(52) 



The scale parameter M is the same as in eq. ([!]), while the coefficient of the 8 function 
is uniquely fixed by consistency with the previously established formulae. We multiply 
eq.(|52|) by x^ and take the Fourier transform. After integrating the derivative by parts 
on the right-hand side, we get the integrals of the known types, eqs.@ and ([3D, while 
the left-hand side becomes just I n , eq.(|J). Thus, we can easily find the necessary 
coefficient. 



The Feynman integrals in the extended differential renormalization, to be compared 



with the expressions in square brackets in eqs.(|46|)-(|50D, can now be evaluated to 
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fig§i 

fig-ie 
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ln 2 (/i 2 /^ 2 ) , 
-InV/z, 2 ), 

hiV/z/ 2 ) -2\n{p 2 /u 2 ) - ln 2 (/i 2 /^ 2 ) + 
(p^Pu/p 2 ) [hiV/z/ 2 ) - 21n(p 2 /z/ 2 ) + 2] , 

InV/z/ 2 ) - 21n(p 2 /^ 2 ) - ln 2 (/i 2 /z/ 2 ) + 2 
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4 9\M> 



- {Pv-PvIp 2 ) hi(p 2 /z/ 2 



(53) 
(54) 
(55) 

(56) 
(57) 



We see that eqs. (|56|) and (|57|) differ from the finite parts of eqs.(fi^) and (|50|) in dimen- 
sional renormalization by the lack of the lower-degree logarithm of /i 2 . As a result, we 
get quite a different charge renormalization 



^diff 
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In 



2) ( 21n 7) ^ + {n ~ 2) (" 3 ln2 7 + 21n 7 - 2 ^ 



h 2 



(4tt 



>2 ! 



/WO = -2(n - 2) /i 2 /(4tt) + • /i 3 . (58) 

Equations ( |5*8D and ( |5*TD disagree beyond the range allowed by the scheme arbitrariness 
in perturbation theory because in a one-charge model the two-loop (3 function must be 
scheme-independent fill . 

Thus, the infrared extension of differential renormalization completely fails in two- 
dimensional models with intrinsic infrared divergencies. The results of multiloop 
renormalizat ion-group calculations may prove to be incorrect or ill-defined. In higher 
orders of perturbation theory, ambiguities cannot be avoided even if we use the same 
scale fi=u for renormalizing infrared and ultraviolet divergencies. 



3 D = 4: Artificial infrared divergences 

There are no off-shell infrared divergencies in D=A. However, they may appear 
as a result of the infrared rearrangement in Feynman diagrams, aimed to the simpli- 
fication of multiloop calculations |J, without which some complicated diagrams may 
technically occur unattainable for analytic computation. Indeed, the usual trick is 
to put some (or all) external momenta to be zero (that is to integrate over the cor- 
responding co-ordinates in x space), and then to calculate divergencies, for example, 
the poles in dimensional regularization. According to the theorem proven in ref. 0, 
after subtraction of sub divergencies, the singular part of a diagram is a polynomial in 
external momenta and internal masses, which is reduced to a constant independent of 
them in the case of a logarithmic divergency. This allows one to simplify calculations, 
taking care of the infrared divergencies that may appear as a result of nullifying some 
momenta. These artificial infrared divergencies can be subtracted by means of the 
infrared R operation 0. 



Our aim here is to check whether the formalism of differential renormalization can 
be extended to perform the infrared R operation in the same way as it performs the 
usual ultraviolet R operation in a renormalizable model. 

To be concrete, consider the massless scalar field theory <P i {D _ i) • The basic equation 



of differential renormalization jl] is the four- dimensional analog of eq. ([!]), 



and its counterpart in momentum space is 



i/(/r 



ln(x 2 M 2 )/x 5 



\n(p 2 /v 2 )/p 2 



(59) 



(60) 



Here, M and v are the ultraviolet and infrared scales, respectively 

Following the rules of the differential renormalization method [Ij], we perform the 
calculations, replacing all the singular expressions according to eqs.(|59D and flBOD , and 
then integrating the derivatives by parts and ignoring surface terms. However, we 
should point out an important difference from our line of reasoning with the set of 
Fourier-transformation formulae in two dimensions. If we allow us to differentiate 
Feynman integrals, generating singular expressions, then the circle of checking con- 
sistency of the integrals does not close by a nonzero additive constant. Thus, we are 
forced to follow a more cautious way, avoiding any explicit differentiation, the result 
of which may be singular. Neither shall we attribute any meaning to the logarithm at 
its singular zero point. 

Obeying these rules and using an intermediate analytic regularization when needed, 
we get the following set of formulae: 
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(61) 
(62) 
(63) 
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where F(a) is the same as in eq.(|]) above. In particular, to the lowest orders, one gets 



4tt 2 /p 2 
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(66) 
(67) 
(68) 
(69) 

(70) 



and the same for the space-time integrals i n and j n with the interchange x 2 <-^p 2 , 
M 2 <r^l/v 2 , fi 2 ^-l/N 2 . The relation between \x and M remains the same as in two 



dimensions; capitals refer to x space. Using eqs.(|5^)-(|64]), d65|)-([?0|), we can perform 
all the calculations. 

In principle, the set of formulae can be supplemented with renormalization-scale- 
independent singular equations 



n x (l/x 2 ) = -16 7T 2 6 4 (x) , 



n p (i/p 2 ) = -16 7T 2 5 4 (p) • 



However, integrating a derivative by parts, we can avoid using them explicitly, for the 
sake of safety always keeping ourselves at least one step off the singular threshold. 

Consider, for pedagogical purposes, the vertex function in the two-loop approxima- 
tion (see fig.||). 
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Figure 6: The vertex function of the 4 theory in the two- loop approximation. 

We demonstrate the infrared peculiarities of the extension of differential renormal- 
ization by the calculation of the renormalization-group (3 function. It coincides with the 
anomalous dimension of four-point vertex up to the two-loop propagator contribution 
ignored hereafter. 



3.1 Non- exceptional momenta 

To simplify the calculation, we reduce the diagrams to the propagator type by 
nullifying some external momenta. One can do this in two different ways. The first way 
is to take non-exceptional external momenta without creating any infrared divergencies. 
Then, according to eqs.(|59| ) - (|69|) , we have [Hereafter, we ignore obvious factors of 7r 2 
coming from the integrals, together with (27r)~ 4 from the loops, so that the contribution 
of each diagram should be divided by (167r 2 ) n , where n is the number of loops.] 
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x a; 

Summing up all these contributions, we get the invariant charge as follows: 

h diS = h+\h 2 ln(p 2 //i 2 ) + h 3 [l lnV/^ 2 ) - 3 ln(p 2 //i 2 ) + 3 
The (3 function is defined by 

2 dh 
(3(h) = fi ^ , 



(71) 



(72) 



while the invariant charge h is /i-independent. Differentiating eq.([fl]) with respect to 



fj, , we get 
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Compare this calculation with the one that uses dimensional regularization and the 
MS scheme. One has 
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Thus, for the bare and invariant charges, we have, respectively, 

h B =(^)'[h+ | (h 2 -h*) /e+ \h?/e 2 



and 



h^ = h+lh 2 Hp 2 /n 
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(74) 



(75) 



We can now calculate the (3 function, differentiating either eq.(|74]) or eq.(75) with 
respect to /i 2 and taking into account eq. ([72]) . Both ways lead to eq.(J73) for the (3 
function. 



3.2 Exceptional momenta in differential renormalization 

The second possibility is to take some exceptional momenta. In our case, setting 
two of them equal to zero, we get the following decomposition (see fig.0). 




Figure 7: Two-loop vertex diagrams in case of exceptional momenta. 

Now, we ought to consistently remove infrared divergencies by defining the artificial 
tadpole graphs which are both ultraviolet- and infrared-divergent. We can use the 
invariance of the regularization, to get a definition, 



d 4 x def 
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where we have evaluated the integrals, according to eqs. (|59|)— (|7Q|) . 



Consider now the two-loop diagrams. We have 
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= - \ lnV//i 2 ) + Hp 2 /ii 2 ) ln(z/ 2 //i 2 ) + \n(p 2 /fi 2 ) - 1 . 

As for the last diagram of fig .[7], it can be calculated in two different ways: first, 
integrating over y in a usual way, and then considering the x tadpole-type integral; or 
conversely. Proceeding in both ways, we get, respectively, 
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Coincidence of these two expressions gives us, thus, a self-consistency check of the 
definition of the artificial tadpole graph (|76|) . 

Combining these two-loop contributions with those already found for non-excep- 
tional momenta, we finally get the invariant charge 
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We see that eq.(|78D differs from eq. (|7T|) . However, differentiating it with respect to 
fi 2 , we obtain the j3 function in agreement with eq. 



3.3 Exceptional momenta in dimensional renormalization 



To get a better understanding of the situation with exceptional momenta, we again 
compare the extended differential renormalization with dimensional regularization in 



the MS scheme. In the latter case, we use the infrared R operation to remove infrared 
divergencies: 
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if.R'i? = o + — 

2s 2 2s 



iT = iln 2 (/i 2 /z/ 2 ) + ln(/i 2 /z/ 2 ) . 



Combining everything together, we get for the bare charge precisely eq.([F4]), and 
for the invariant charge 



UZ =h + h 2 [ln(p7// 2 ) + \ ln(z/ 2 //i 2 ) - 2] + h 3 [ § In 2 (p 2 //* 2 ) 
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(79) 



Note that eq . (|79|) does not coincide with eq. (|78|) . However, if we differentiate eq.([j 
with respect to /i 2 , we reproduce the correct value of the (3 function eq. (|73|) . 



Hence, in dimensional regularization, when one uses R* operation, the f3 function 
calculated both from infinite and finite parts is independent of external momenta in 
the MS scheme. This is also true for the finite expressions obtained with the aid of the 
extended differential renormalization, which give us the correct value of the (3 function 
despite the fact that the finite corrections to the invariant charge eqs.(frT|) and ( f78|) are 
different from those in the MS scheme, eqs. (|75|) and (|79|). This difference is, however, 
trivial and reflects just the renormalization scheme arbitrariness. Moreover, to two 
loops, the finite part of each diagram in the extended differential renormalization can 
be obtained, up to a constant, from the corresponding MS renormalized expression by 
a shift: ln/z 2 — dn/x 2 — 2. This statement happens to be independent of the presence or 
absence of infrared divergencies. 



3.4 Higher-order ambiguities 

This pleasant picture is, however, spoiled in higher orders. To see this, let us 
consider the following four-loop artificial tadpole graph: 

d A xd 4 yd 4 zd 4 t 




(y 2 ) (y - x) 2 (x - zf [(z - t) 



i 2 



To evaluate the internal propagator-type subgraph, we need to add a new equation to 
the set of basic formulae, eqs.fl59l)-(fr0|), 



1 



ln(a; 2 M^ 



n* 

32 3 



x 



+ — 7T 2 n x 6i(x) , 
lb 



10) 



where the coefficient of the last term is strictly fixed by the requirement of consistency 
with the rest of the formulae after multiplying both sides by x 2 . Performing the internal 
integration with the aid of eq. fl8"0|) , we reduce the diagram to the following integral, 



d 4 x d 4 y [ln(x 2 M 2 ) + § 



(y 2 ) (y - x) 2 x 2 



51) 



which, again, can be integrated first over y and then over x, or in the opposite order. 
We have, respectively, 



dV 



\n(x 2 M 2 ) \n(x 2 M 2 ) + 



d 4 x 
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ln(x 2 M 2 ) [ln(x 2 M 2 ) + | 



x 2 — 2xz + z 2 



(x - z) A 



I ln 3 (^ 2 //i 2 ) + \ ln 2 (^ 2 /^) - f m^V/i 2 ) + | C(3) , 



(82) 



d 4 g . ln(g 2 //, 2 



d'y 
(/J 

\ I 4h { i ln2 ^ 2]v2 ) + [i - ln (^ 2 )] Hv 2 n 2 ) + 1} 



InV/V) + I InV/V) - f ln^ 2 /^ 2 ) " |C(3) - 1 



In contrast to the previous case of eqs.([77D, the two expressions (p2|) and (|82f ) are 
different. It turns out to be a general situation, which simply has not yet manifested 
itself at the two-loop level but becomes evident at the fourth loop [In fact, integrals 
like eq.QSip appear already at the three-loop level]. Thus, we see that, in D=4 either, 
the infrared extension of differential renormalization does not provide us with a self- 
consistent definition of tadpole graphs. On the other hand, abandoning the infrared R 
operation and just nullifying artificial tadpoles, we would not be able to make use of 
infrared rearrangement at all. 

Therefore, we conclude that the infrared extension of differential renormalization 
cannot perform the infrared R operation, like the original scheme does the ultraviolet 
R operation. 



4 Summary 

We have examined the possibility of generalizing differential regularization in an 
invariant fashion to theories with infrared divergencies. Both in D=2 and D=4, the 
basic differential identities of the method, written in co-ordinate and momentum space, 
lead to a definite set of consistent formulae for divergent Fourier integrals regularized 
by two scale parameters, /1 and u, which remove ultraviolet and infrared singularities, 
respectively. 

The principle of invariance of the regularization allows us also to derive definitions 
for the tadpole-type integrals (without external momenta or co-ordinates), which in- 
termix infrared and ultraviolet divergencies, by splitting the integrals up into a sum 
of separate infrared, ultraviolet, and regular items. However, the values for tadpole 
graphs of a complicated structure in higher orders of perturbation theory prove to be 
ambiguous, depending on the order of evaluation of their subgraphs. This ambiguity 
cannot be fully eliminated even if we try to somehow relate v with /j, (which makes 
low-order tadpole diagrams to be zero, but nevertheless, does not lead to a unique 
determination of higher-order graphs; moreover, any information about ultraviolet log- 
arithms is completely lost then). 

In two-dimensional a models, where intrinsic infrared divergencies at zero masses 
are present, the results of the extended differential renormalization disagree with the 
finite minimally renormalized results of dimensional regularization. The difference goes 
beyond the range allowed by renormalization-scheme arbitrariness in perturbation the- 
ory. For the 0(n) sigma model with a single coupling constant, the two-loop coefficient 
of the j3 function turns out to be zero when we use the infrared extension of differen- 



tial renormalization. This contradicts the old well-established result which should be 
scheme-independent. 

As concerns the use of the infrared R operation in dimensional regularization, we 
have verified by direct two-loop calculations that in renormalizable theories (n-field in 
D=2—2e and 4> 4 in D=4—2e) it successfully recovers ultraviolet logarithms in the finite 
parts of the diagrams with exceptional momenta, so that the (3 function, calculated 
through the finite invariant charge, proves to be independent of any logarithms (of 
momenta and renormalization scales). It coincides with the result derived in the con- 
ventional way from the coefficients of the first singularity in e, which does not depend 
on the choice of the momenta for logarithmically divergent diagrams. 

The two-loop results of the infrared extension of differential renormalization in 
the </>(£> =4 ) theory agree with dimensional renormalization up to finite counterterms. 
However, higher-order logarithmically divergent artificial tadpole diagrams (resulting 
from infrared rearrangement) suffer from an irremovable ambiguity. Thus, calculations 
with exceptional momenta cannot be performed consistently above two loops. 

Our final conclusion is that the program of constructing an invariant generalization 
of differential regularization and renormalization, to deal with infrared divergencies, 
has failed. 
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